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'■ Abstract 



Three dimensional unsteady flow of fluids in the Lagrangian description is 
considered as an autonomous dynamical system in four dimensions. The 
condition for the existence of a symplectic structure on the extended space 
is the frozen field equations of the Eulerian description of motion. Integral 
invariants of symplectic flow are related to conservation laws of the dynamical 
equation. A scheme generating infinite families of symmetries and invariants 
is presented. For the Euler equations these invariants are shown to have a 
geometric origin in the description of flow as geodesic motion; they are also 
interpreted in connection with the particle relabelling symmetry. 
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1 Introduction 



In spite of the powerful analytical tools in the study of fluid dynamics, it has 
been argued that the three-dimensional Euler equation requires more than 
just better functional analysis. The need for some geometry, more specifically, 
the necessity of symplectic geometry for tackling some open problems of Euler 
flow is stressed in 

Symplectic techniques have been used, on one side, to cast the Eulerian 
dynamical equations into Hamiltonian (Lie-Poisson) form@ and, on the other 
side, the Lagrangian description of steady fluid motion in two and four di- 
mensions has been shown to have the structure of a completely (Liouville) 
integrable Hamiltonian systemS. 

In this work we shall be concerned with the kinematical problem of un- 
steady fluid motion, namely, the Hamiltonian structure of the velocity field 
in the Lagrangian description of motion. This will be an application in co- 
variant framework of the construction in @] of the Hamiltonian structure of 
time-dependent dynamical systems. The characterization of unsteady three 
dimensional flow of a fluid as a Hamiltonian system in the symplectic frame- 
work will involve the dynamical equations for frozen-in fields, thereby mani- 
festing an interplay with the Eulerian description. 

The covariant symplectic techniques will enable us to construct Lagrangian 
and Eulerian invariants as well as to elucidate the connection between them. 
We shall show that an infinite family of invariants can be constructed for 
any hydrodynamic system admitting a single vectorial frozen-in field. For 
the Euler equation the geometric origin of this family will be found in the 
description of flow as a geodesic motion on an infinite dimensional group. 
We shall argue that the associated infinitesimal generators belong to the left 
Lie algebra of particle relabelling symmetry. 

2 Motion of Fluids in Three Dimensions 
2.1 Kinematic Representation 

Let the open set D C R 3 be the domain occupied initially by the fluid and 
x(t = 0) = Xq G D be the initial position, i.e., a Lagrangian label. For 
a fixed initial position xq, the Eulerian coordinates x(t) = gt{xo) define a 
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smooth curve in R 3 describing the evolution of fluid particles. Thus, for each 
time t G R, the orientation (or volume) preserving embedding g t : D — > R 3 
describes a configuration of fluid. A flow is then a curve 1 1— > g t in the space 
of all such maps. The time-dependent Eulerian (spatial) velocity field u t that 
generates g t is defined by 

^ = = ( u t °9t)(xo) = u(t,x) (1) 

where Ut o ^ is the corresponding Lagrangian (material) velocity field@ . 

Since g t is volume preserving, u t (x) is a divergence- free vector field over 
R 3 and eq(]B) is a non-autonomous dynamical system associated with it. 
Equivalently, ([!]) can be represented as an autonomous system by the vector 
field 

d t +u{t,x), m = u-V (2) 

on the extended space R x R 3 . 

The Lagrangian description of fluid motion is the study of the system 
of non-autonomous ordinary differential equations (]!]), or, equivalently, its 
autonomous counterpart associated with fl2|). We shall be interested in its 
Hamiltonian formalism which will turn out to imply and be implied by dy- 
namical equations for frozen-in Eulerian fields. 



2.2 Dynamical Equations 

The dynamics of the the Eulerian velocity field u(t, x) is determined by 
Newton's law. For an ideal incompressible isentropic fluid in a bounded 
domain D with normal n this is given by the Euler equations 

^ + U . VU ^-V(IP) (3) 

where p is the density and P is the pressure. The additional conditions 
V ■ u = and u • n|ao = on the velocity field characterize it as the 
generator of volume preserving diffeomorphisms. The identity u ■ Vu = 
|V|u| 2 — u x (V x u) can be used to bring the equation (|3]) into the form 

^-ux(Vxu) = -V(-P + ^) (4) 
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which we shall refer to as the rotational Euler equation. In terms of the 
vorticity vector w = V x u the Euler equation becomes 

— - V x (u x w) = (5) 

and for our purposes we shall assume the existence of a scalar field h advected 
by the fluid 

^ + u-Vh = (6) 

which further satisfies w • V/i 7^ 0. This latter condition does not affect the 
generic structure but is only necessary to present the results in covariant 
symplectic form. 



3 Symplectic Structure of Kinematic Descrip- 
tion 

A symplectic structureB'H' on a manifold N of even dimension 2n is defined 
by a closed, non-degenerate two-form Q. It is exact if there exist a one-form 
9 such that Q = —d6. Darboux's theorem quarantees the existence of local 
coordinates {q l ,Pi) i = 1, ...,n in which Q has the canonical form dq l A dpi. 
The 2n— form (— l) n Q n /n\ is called the Liouville volume. A vector field X 
on N is called Hamiltonian if there exists a function h on N such that 

i(X)(n) = dh (7) 

where i(X)(-) denotes the inner product with X. The identity i(X)(dh) = 
which follows from (|7|) is the expression for conservation of h under the flow 
of X. With the correspondence ([!]) between functions and vector fields, the 
Poisson bracket of functions on iV defined by 

{f,g} = Q(X f ,X 9 ) = n-\df,dg) (8) 

satisfies the conditions of bilinearity, skew-symmetry, the Jacobi identity and 
the Leibniz rule. This enables us to write the dynamical system associated 
with the vector field Xh in the form of Hamilton's equations 

^ = {*,/*}. (9) 
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Proposition 1 (1) The suspended vector field dt + u on Rx R 3 is a Hamil- 
tonian vector field with the symplectic two-form 

fi e = -(Vh + u x w) ■ dx A dt + w ■ (dx A dx) (10) 

and the Hamiltonian function h. 

(2) The potential vorticity — w- Vh = ph is the invariant Liouville volume 
density of symplectic structure. 

(3) fl e is exact with the one-form 

-6 e = if) e dt + u-dx, -i/j e = h + -P + -u 2 . (11) 

P 2 

(4) In Darboux coordinates, Q e = dq A dp + dt A dh can and the velocity 
field has the Clebsch representation 

u = Vs + pVq (12) 

where s is the generating function of the canonical transformation. 

Proof: Q e is closed via eq([|) and the divergence free property of vorticity. 
The non-degeneracy is a consequence of the assumption that — w ■ Vh = ph ^ 
where ^ 

~Q e A fi e = p h dx A dy A dz A dt (13) 

is the definition of the Liouville volume. Its invariance can either be checked 
directly by showing that 

div Ph (d t + u) = ^ + V- (p h u) = . (14) 

or more geometrically, follows from the recognition that eqflU), when written 
in terms of differential operators u = u • V and ui = w • V as 

^+[u,uj] = -(V-uV , (15) 

is the criterion for u to be a time- dependent symmetry of u. Then, the fact 
that, for dynamical systems, action of symmetries on conserved quantities 
gives conserved quantities@ implies the invariance of ph- See also |§ for a 
physical description of Lagrangian invariants produced in this way. 
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Eq(^) is used to write dt + u in Hamiltonian form. The definition of 
vorticity field as a curl and the rotational Euler equation (f|) turns Q e into 
an exact symplectic form with (|TT|). Clebsch variables for velocity field as 
well as the form of Hamiltonian function 

Os do 

h mn = ^ e --- P - (16) 
in Darboux coordinates follows easily from the defining relations 

ds = 6 can - 9 , -6 C an = pdq + h can dt (17) 

for these coordinates. • 

An important conclusion of this result is that the Eulerian dynamical 
equations determining the form of velocity field enforce its suspension to 
be a Hamiltonian vector field. In the next section we shall show that with a 
slight modification the physical framework of proposition (1) can be extended 
to other hydrodynamic systems. 



4 Generalizations on Dynamical Equations 

The conditions for the existence of Hamiltonian structure of the velocity field 
can be made independent of the dynamics of the fluid itself by replacing vor- 
ticity in Q e with a divergence free field which is advected by fluid motion. 
This allows one to consider fluid motion under more general forces. Velocity 
fields governed by the Euler equation (|3]) has already such a generalization 
which involves little alteration of their Hamiltonian structure. Namely, any 
lift force can be incorparated to the right hand side of eq(|3|) by the replace- 
ment V/i + u x w i— > Vh + u x w + F in Q e . The foremost example of this 
sort of force fields with vanishing component along u are the magnetic forces 
acting on electrically conducting fluids. 

Proposition 2 Let the dynamics of u be governed by 

^ + u-Vu = F (18) 
and assume that the divergence-free field B and the function if satisfy 

^-Vx(uxB)=fl, ^ + U . V ^ = (19) 
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which are the frozen field equations. Then 

(1) dt + u is a Hamiltonian vector field with the symplectic two-form 

n = -(Vy? + u x B) • dx A dt + B ■ (dx A dx) (20) 

and the Hamiltonian function if. 

(2) ptp = — B ■ Vip is the invariant Liouville volume density. 

(3) If moreover B = V x A for some vector potential A, then Q is exact 

tt = -d6, -6 = ipdt + A-dx (21) 

where if) is determined by the equation 

dA 

--ux(VxA) = V(^^). (22) 
at 

(4) The Darboux coordinates for Q defines the Clebsch representation of 
the vector potential A. 

We shall now show with an example from dynamo theory that the Hamil- 
tonian character of Lagrangian description is determined by dynamics of the 
B field rather than the velocity itself. The dynamics of the velocity field 
is determined by the so-called magneto-convection equation fll8| ) where the 
complicated velocity dependent force field per unit mass 

F = ~B x (V x B) + zA7 2 u-2fi x u + Cg - vru (23) 

includes Lorentz, viscous, coriolis, gravity forces and the reduced pressure. 
We refer to for a description of the physical framework and terminology 
involved. The divergence free magnetic field satisfies 

— - V x (u x B) = r?V 2 B (24) 

where rj is the magnetic diffusivity. With this equation the two-form 

tt d = -(Vp + u x B - r]V x B) ■ rfx A dt + B ■ (dx A dx) (25) 

is symplectic. It is exact with the one-form fl2"T| ) where the scalar potential 
ip is determined from eq(^) whose right hand side is replaced by V(ip + 
ip — f]V ■ A) + f]V 2 A. However, the existence of a Hamiltonian function 
is prevented by the diffusion term in eq(ES). For rj = 0, or equivalently for 
infinite conductivity, the suspension of velocity field is Hamiltonian regardless 
of how complicated the force field may be. 
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5 Integral Invariants and Eulerian Conserva- 
tion Laws 



Conservation laws of dynamical equations are divergence expressions of the 
form 

dT 

— + V • P = (26) 

where the conserved density T and the flux P are functions of t, x, the Eu- 
lerian fields and their derivatives. Some of these Eulerian conservation laws 
are known to be related to the Lagrangian conserved quantities. In this sec- 
tion we shall use the symplectic structure to construct invariant differential 
forms of the Lagrangian picture and show with examples their connection 
with conservation laws of the Eulerian description. 

We start with the observation that eq(^) can be expressed as the closure 
of a three-form in four- dimensional space-time. If is the closed three-form 
corresponding to eq(^), then we have for its Lie derivative 

£ dt+u {6) = dS (27) 

where S = i(d t + u)(Q). This means that G is a relative integral invariant of 
dt + u. Conversely, it turns out that some conserved densities of the Eulerian 
dynamics can be constructed out of invariant forms of symplectic structure 
of Lagrangian kinematics. The formulation of the continuity equation (|14]) 
as the closure of invariant three- form 0a/ = pl\^=i(dx t — u l dt) can be viewed 
as a primitive example of the sort of relations we wish to consider. 

5.1 Helicity Conservation 

The canonical one-form 9 is a relative integral invariant for its Lie derivative 
£g t+u (9) = d(i(d t + u){9) — (p) is only closed. Its exterior derivative which 
is the symplectic two-form is then an absolute invariant. Using them we 
construct the (relative) invariant three-form 

Q H = 9 A d6 = A • V x A dx A dy A dz + 

((V> + u ■ A)B - (A • B)u - A x W) ■ (dx A dx) A dt (28) 

which consists of the well-known (magnetic) helicity and corresponding flux 
densities. The identity d(8 A d6) — Q A Q = expressing the fact that its 
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closure is an absolute invariant, gives the conservation law 

-^(A • B) + V • ((A • B)u - (V> + ip + u ■ A)B) = (29) 

for the dynamics governed by general force fields. In the case of Euler flow 
eq (|29|) is the conservation law for the helicity u • w which is the Casimir of 
the Lie-Poisson bracket for the Euler equations in the vorticity variableS. 

5.2 Energy Conservation 

The energy conservation arises from a somewhat different consideration. If we 
regard the closure conditions of symplectic two-form Q as the homogeneous 
part of Maxwell equations, the inhomogeneous part will define the exact 
three-form 

® E = q dx A dy A dz - J ■ (dx A dx.) A dt (30) 
where the expressions for charge and current densities 

V-E = g, -— + VxB = J (31) 
at 

are the usual form of inhomogeneous Maxwell equations. From Q we read 
off the electric field E = — (Vip + u x B). The closure of 9^ is equivalent to 
the well-known Eulerian conservation law 

da 

de E = (-^ + V ■ J) dt Adx A dy Adz = (32) 

which easily follows from definitions (|3~TD. For the symplectic structure ( p0|) 
the conserved density q is given by 

q = -VV-w-B + u- (V x B) (33) 

which reduces to V 2 ip in the Coulomb (or transverse) gauge V • A = 0. Note 
that this gauge condition is satisfied automatically for 9 e . In the case of Euler 
equation we have 

q e = _V 2 /i -w 2 + u-(Vxw) = -V 2 (/i + -P+ -u 2 ) = V 2 ^e (34) 

P 2 

which is the form of energy density analogous to the one considered in || for 
the case of steady flows. 
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6 An Infinite Family of Invariants 



The introduction of is more crucial for the present section than for the 
relatively trivial conservation law (|32|) . It naturally exists via definition (j3ll) 
for any hydrodynamical system having one vectorial frozen-in field. One of 
the obvious exceptions is the viscous flow of the Navier-Stokes equation, and 
another is the fluid with finite conductivity which we have discussed before. 

6.1 A New Symmetry and Associated Invariants 

We shall show that the requirement on Qg to be an absolute invariant results 
in a new symmetry of the velocity field, gives an Eulerian conservation law 
and turns q into a conserved density of the Lagrangian description as well. 

Proposition 3 Qe is an absolute invariant of dt + u if and only if the 
divergence-free vector field qu — J is a time- dependent symmetry of u. In 
this case, 

(1) q is also a Lagrangian invariant and 

(2) A - (qu — J) is an Eulerian conserved density depending on dynamics 
of velocity field through the definition of 3 . 

Proof: Since the Lie derivative of Qe 

£ dt+u (6 E ) = d((qu - J) ■ (dx A dx) - (J x u) ■ dx A dt) (35) 

is exact, it is in general a relative invariant. The condition that Gg is an 
absolute invariant is the vanishing of the right hand side of ([35]) . This implies 
that qu — J must be divergence-free and satisfy 

9{qU m J) + [u, qu-J] = -(V • u){qu - J) (36) 

which is the equation defining a symmetry of u. In particular, given u and 
q for an incompressible flow, the solutions J of linear inhomogeneous first 
order equations 

dJ r Tl du 
makes qu — J into a symmetry of u. 
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The invariance of q under the unsteady flow of u 



-g + u ■ Vq = (38) 

follows from the divergence-free condition and the conservation law (j32|). 

To prove (2), we first observe that if Qe is an absolute invariant, so is 
the closed two-form = i(d t + u)(Qe)- Then the three-form 9 A is a 
relative invariant and the conservation law 

J^(A • (qu - J)) + V • (A x (J x u) - {ip + ^)(qu - J)) = (39) 

follows from the identity d(9 AS^)+f2AS£; = 0. • 



6.2 New Symmetries and Underlying Geometry 

It can easily be shown that the Lie bracket of symmetries B (or w) and 
qu — J of u is a new time-dependent symmetry. In fact, this is true for 
any two symmetries of u. Thus, repeated applications of the Lie bracket 
generates an infinite dimensional algebra of time-dependent symmetries of 
the velocity field. 

New integral invariants can then be constructed from Lie derivatives and 
inner products of known differential invariants with these symmetries (see for 
example, exercise 6.18 of ref ||). Relations with Eulerian conservation laws 
are established by taking appropriate products and combinations of invariant 
forms. The derivations of conservation laws (|29"D and (|39"D are examples of 
this construction. 

To this end, we wish to obtain a geometric characterization of the origin of 
conservation laws generated by the above procedure. In order to do this, we 
shall first restrict them to a particular gauge in which the densities are also 
conserved along the trajectories of velocity field. The resulting Lagrangian 
conservation laws will be familiar from geometry. The appropriate gauge can 
be characterized by demanding the canonical one-form 9 to be an absolute 
invariant of dt+u, that is, by the condition ip+i/j + u-A = constant. Then for 
an incompressible fluid, the conserved densities in eqs(p9|) and (|39f), namely 
the magnetic helicity A • B, and A • (qu — J) are also Lagrangian invariants. 

We shall now show by further restricting to Euler equation that the in- 
finite family of invariants have a geometric (rather than dynamic) origin in 
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the description of Euler flow as a geodesic motion on an infinite dimensional 
group. For the Euler equation (H) the gauge condition implies that 6 e is 
absolutely invariant on surfaces defined by u 2 /2 — P(p)/p = constant. The 
corresponding Lagrangian invariants are the helicity u • w, and u • (gu — J). 
The restriction on dynamics for which the velocity field admits this sort of 
invariants is contained in the following straightforward result. 

Proposition 4 Let v be a time- dependent symmetry of velocity field u whose 
dynamics is governed by eq ^H{ ). If 

v • (F + ^Vw 2 ) = (40) 
then u ■ v is a Lagrangian invariant. 

This explains the relevance of the chosen gauge for the Euler equation. As an 
immediate consequence, there follows the conservation of u - B in a magnetic 
relaxation problem which was also indicated in ||. Note that there is no 
dynamical restriction for invariants of the form A ■ v. 

The Lagrangian conservation laws which are in the form of a product of 
velocity with a symmetry are familiar from geometry where one proves, as 
for example in |TD[ , that if u is tangent to a geodesic curve on a Riemannian 



manifold and if v is a Killing vector (an infinitesimal isometry), then u ■ v 
is an invariant of any motion with velocity u. The analogy with the above 
invariants of Lagrangian description will become transparent if we recall that 
in Lagrangian coordinates the Euler flow is a geodesic motion on the infinite 
dimensional group of volume preserving diffeomorphismsQ'@. 



Since Ut is right invariant (c.f. eq(|l])), the motion is generated by left 
translation. The symmetry criterion for v t (x) = v(t,x) is the infinitesimal 
expression for the invariance of v t under the pull-back g*(v t ) = Tg^ovtogg by 
geodesic motion. So, any symmetry generator is a left invariant vector field. 
Moreover, since the starting symmetries w and qu — J are divergence-free, it 
follows from the identity 

div ([v , w)) = v div{w) — w div{y) (41) 

that they all are divergence-free and hence are elements of the algebra of 
volume preserving diffeomorhisms. The left invariance of symmetries implies 
that they are the generators of the right action or, in fluid mechanical context, 
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the particle relabelling symmetry. Thus, the conserved densities we have been 
considered so far are all associated with this general symmetry property of the 
Euler flow. It has indeed been shown explicitly (see |T2|] and the references 
therein) that many of the well-known conserved densities of fluid motion are 
derivable from a so-called general vorticity conservation law connected with 
the particle relabelling symmetry. We may therefore claim that the infinite 
dimensional symmetry algebra constructed above is the left Lie algebra of 
group of volume preserving diffeomorphisms and it consists of infinitesimal 
generators of the particle relabelling symmetry. 



7 Discussion and Conclusions 

We have shown that the symplectic structure of kinematic description presents 
an interplay with dynamical description. The connection between invariants 
of velocity field and conservation laws of dynamical equations is another 
manifestation of this interplay. 

We obtained the Eulerian dynamical equations (H), ( |19D and fl22|) from 
a purely geometric treatment inherent in the Lagrangian description. In an 
earlier work, described in [131, Kozlov obtained eqs(f|) and (H) in constructing 



integral manifolds of finite dimensional non-autonomous dynamical systems 
thereby first recognizing the connection with hydrodynamics of an ideal fluid. 
Motivated by physical arguments Kuz'min introduced in [[| (see also ||) 
fields satisfying (in a particular gauge) eq(p2|) and used them to construct 
Lagrangian invariants. 

In this work we presented a much more suitable construction relying on 
symplectic geometry for a systematic and exhaustive investigation of both 
Lagrangian and Eulerian invariants. Moreover, the construction inherits geo- 
metric settings not only for the invariants themselves but also for their defin- 
ing equations. We therefore conclude that the Lagrangian description, which 
is seldom used in fluid dynamics, provided with a geometrical framework, 
may become a powerful tool in analysis and the construction of invariants of 
hydrodynamical models. 
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